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Qh| Abstract 

When the neutrino mass eigenstate vi is trimaximally mixed, the mixing matrix 
is called trimaximal. The middle column of trimaximal mixing matrix is identical to 
tri-bimaximal mixing and the other two columns are subject to unitarity constraints. 
This corresponds to a mixing matrix with four independent parameters in the most 
general case. Apart from the two Majorana phases, the mixing matrix has only one free 
parameter in the CP conserving limit. Trimaximality results into interesting interplay 
between mixing angles and CP violation. A notion of maximal CP violation naturally 
emerges here: CP violation is maximal for maximal 2-3 mixing. Similarly, there is a 
natural constraint on the deviation from maximal 2-3 mixing which takes its maximal 
<^ ! value in the CP conserving limit. 



X 1 Introduction 

The historical discovery of neutrino oscillations constitutes the only indication of physics 
beyond the standard model. In the absence of a basic understanding of fermion masses and 
mixings, it is natural to seek simplified patterns in the mixing matrix which are consistent 
with the present experimental knowledge and provide guidance for the future experimental 
searches. Such phenomenological patterns are, at best, first approximations to the experi- 
mental data and the next step is to study deviations from these patterns. 

A successful phenomenological ansatz for neutrino mixing matrix which is consistent with 
the present neutrino data [1] was proposed by Harrison, Perkins and Scott [2] and is given 
by 
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In this mixing scheme, the second neutrino mass eigenstates V2 has trimaximal character as 
it arises from maximal mixing of the three flavor eigenstates (z/ e , v^ and u T ). Moreover, the 



third mass eigenstate v% has bimaximal character as it results from maximal mixing of two 
flavor eigenstates viz. v^ and u T . Hence, this mixing scheme is called tri-bimaximal (TBM) 
mixing and gives vanishing #13 and maximal 823 . The mixing angle #13 vanishes in TBM 
mixing because of bimaximal character of z/3, a feature needs to be abandoned in order to 
allow for a non-vanishing #13 and, hence, the possibility of CP violation. However, there is 
no need to abandon the trimaximal character of the z/ 2 in this generalization of TBM mixing. 
Such a mixing matrix, called trimaximal (TM) [3j, will be of the form 
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where the first and third columns are subject solely to unitarity constraints and the middle 
column has trimaximal character. The definition of trimaximality adopted here requires 
that magnitudes and phases of all the three elements in the middle column of the mixing 
matrix are equal. The common phase of the middle column can be rephased away by phase- 
redefinition of the lepton fields. 

The relationship between TBM mixing and TM mixing needs also to be understood from the 
standpoint of the symmetry groups associated with these mixing schemes. The TM mixing 
follows from a neutrino mass matrix which can be parametrized as 
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in terms of four complex parameters a, b, c and d having v4 4 symmetry in the flavor basis 
where the charged lepton mass matrix is diagonal [3]. The matrix is called magic mass 
matrix since its rows and columns all add up to a + b + c and hence one of its eigenvectors 

is (l/\/3 1/VS 1/V3J within an overall multiplicative phase [3]. The above mass matrix 
becomes fi-r symmetric as well if the number of free parameters are reduced to three by 
demanding b = c. The mixing matrix become TBM and the symmetry group enlarges to S4 
after making the above fine tuning. Thus, the neutrino mass matrix having A4 symmetry 
and bringing about TM mixing contains one more complex parameter than the mass matrix 
having S4 symmetry and producing TBM mixing which can be used to accommodate a non- 
vanishing 6*i3 and a Dirac-type CP-violating phase 5. This additional parameter needs to 
be fine tuned to recover the TBM mixing in a mass model based upon A4 symmetry. The 
neutrino mixing matrix will be TM in absence of any such fine tuning of parameters pQ. It 
becomes imperative to analyse the phenomenology of TM mixing scheme in this context. 

The CP violating structure in trimaximal mixing is an interesting feature to study because 
of the mathematical restriction that one of the eigenvector of the mass matrix be (1 1 1) 
which severely restricts the elements of the matrix to coordinate in a way so that they 
produce a real eigenvector. This coordination results into a unique correlation between CP 
violation and the mixing angles. 

A notion of maximal CP violation automatically emerges when unitarity constraints are 
combined with trimaximality. The Jarlskog rephasing invariant J [5] takes its maximal 



value -^y= | t/ii 1 1 U13 1 when #23 is maximal. Similarly, the deviation from maximal 2-3 mixing 
is proportional to #13 and is naturally restricted to lie within the present experimental range 
by unitarity for values of #13 within its experimental bounds [1] . This deviation is zero when 
CP violation is maximal and is maximal when there is no CP violation. 



2 Mixing angles and CP violating phases 

Trimaximality, as defined in Eq. (2), equates magnitudes as well as phases of the three 
elements of the middle column of U. Although, the equality of the magnitudes is preserved 
after the phase redefinitions of the lepton fields, same is not true for the equality of the 
phases. Hence, the definition of trimaximal mixing adopted here is not rephasing invariant. 
This is quite natural since a symmetry may prefer a basis, called symmetry basis, for its 
definition. So, it becomes necessary to find the consequences of TM mixing in terms of 
rephasing invariant quantities. The mixing angles and CP violating phases must be deduced 
from these rephasing invariants only. 

There are four independent CP-even quadratic invariants [6] which can conveniently be cho- 
sen as U^Un, U* 3 Ui3, U21U21 and C/I3C/23 and three independent CP-odd quartic invariants 

m 

J = Im(U ll U^U; i U 22 ), (4) 

h = Im[(U* u U 12 ) 2 } (5) 

and 

I 2 = Im[(U* u U l3 ) 2 }. (6) 

The Jarlskog rephasing invariant J [5] is relevant for CP violation in lepton number conserv- 
ing process like neutrino oscillations and I\ and I 2 are relevant for CP violation in lepton 
number violating processes like neutrinoless double beta decay. 

The consequences of unitarity for TM mixing will be worked out in terms of the above 
rephasing invariant quantities viz. \Uu\, | C/13 1 , \U 2 i\, 1 1/23 1 , J, I\ and I 2 . However, it will 
be rewarding to express the phenomenological implications of trimaximality in term of the 
usual parameters 6*12, #23, #13, a, (3 and 5 of a widely used parametrisation [HE] 
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where s^ = sin% and Qj = cos%. This will be achieved through the relations \Un\ = c\ 3 ci 2 , 

\Ul 3 \ = S13, I C/23 1 = C13S23, 

h = 4 3 c 2 12 s 2 l2 sm2a } (8) 

I 2 = c 2 13 c 2 12 s 2 13 sm2/3 (9) 

and 

J = cl 3 c 23 c 12 s 12 s 23 s 13 sin 5. (10) 



However, some of the relationships will be more apparent when written in terms of rephasing 
invariant quantities themselves. 

Like fi — t symmetry and TBM mixing, the TM mixing is mass-independent texture putting 
no constraints on neutrino masses and Majorana phases a and (3 [8]. Hence, trimaximality 
will have no physical consequences for the Majorana phases. 



3 Strong and weak trimaximality 

A naive way to work out the consequences of unitarity for trimaximal mixing would have 
been to compare the magnitudes of the elements of trimaximal mixing matrix of Eq. (2) 
with the unitary mixing matrix of Eq. (7). This would mean that | C/21 1 , It/22 1 and | C/23 1 
in Eq. (7) are all equal to 4= with no restriction on their phases. In this version of TM 
mixing, trimaximality is imposed only in a weak sense as the elements of the middle column 
of U in Eq. (7) have unequal phases for any non zero value of S and #13. This weak form 
of trimaximality has been already analysed in the literature [H] for its phenomenological 
consequences. If a similar analysis to is to be performed for the trimaximality defined by Eq. 
(2), called strong trimaximality to distinguish it from weak trimaximality, which demands 
that magnitudes as well as the phases of C/ 2 i, U22 and U23 are equal, one is required to 
compare the trimaximal U of Eq. (2) with the unitary U of Eq. (7) only after making its 
middle column completely real by making a phase rotation on the lepton fields. This process 
of making a basis transformation on a mixing matrix to the symmetry basis where the TM 
mixing is defined is quite natural if the definition of symmetry is not rephasing invariant and 
prefers some basis. As stated already, the strong trimaximality is not a rephasing invariant 
concept unlike weak trimaximality which is a rephasing invariant notion. However, instead of 
adapting a unitary matrix by phase rotation to match it with trimaximal mixing as outlined 
above, the present work attempts to obtain an unitary parametrisation which is trimaximal 
in the strong sense from the very outset. 

The strong trimaximility implies that the mass eigenstate V2 is a democratic and coherent 
combination of the three flavor eigenstates whereas the week trimaximality relaxes the con- 
dition of coherence allowing different flavor eigenstates to mix with equal amplitudes but 
mismatching phases. However, the distinction between strong and weak trimaximality is 
not merely conceptual, it is also phenomenological in nature having distinct predictions for 
CP violation. The weak trimaximality is less interesting as it doesn't fully utilise the phase 
restrictions on the mass matrix coming from A 4 symmetry. 

A mixing matrix which is trimaximal in the strong sense was first proposed and studied 
phenomenologically by Bjorken, Harrison and Scott (BHS) [10] in which the only free complex 
parameter was U13. A neutrino mass matrix having magic symmetry in the basis of a diagonal 
charged lepton mass matrix was prescribed by Lam [TT] which gives rise to TM mixing in BHS 
parametrization. In another attempt for the minimal modification to the TBM mixing, He 
and Zee [12J obtained a mixing matrix, which is also trimaximal in strong sense, containing 
only one free complex parameter. The neutrino mass matrix in a model proposed by Grimus 
and Lavoura [13] is magic and hence gives rise to TM mixing in the strong sense. However, 



the mass matrix contains three complex parameters and therefore is more restrictive than a 
most general magic mass matrix [Eq. (2)] containing four complex parameters. Similarly, 
the resulting TM mixing matrix also contains lesser number of parameters than the most 
general case. All the above forms of TM mixing studied phenomenologically in the literature 
are more constrained than the most general form given by Eq. (2) and can be expressed as 
the product of the TBM mixing matrix of Eq. (1) by a generalized 1-3 rotation of the type 
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from the right. However, the resulting mixing matrix is not the most general TM mixing 
matrix. 



4 Unitary constraints 



The TM mixing matrix of Eq. (2) can be written as 



U 



( 


U u 


i 

v/3 




u 21 


1 
v/3 


\ 


-U n - U 21 


1 



u l3 

u 23 

-Uis — U 23 



:n) 



by straightforwardly eliminating two elements viz. U 3 i and U 33 using the two orthogonality 
conditions for the three columns. The remaining 4 elements are further constrained by the 
following seven constraints coming from the unitarity condition UlP — U*U — 1: 
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(13) 
(14) 

(15) 

(16) 
(17) 

2Im{U n U* 3 ) + Im{U 2l U* lz ) + Im{U n U* 3 ) + 2Im{U 21 U* 3 ) = 0. (18) 

Eq. (11) is a unitary parametrisation of TM mixing subject to the above seven unitarity 
constraints. All the parametrizations of TM mixing proposed in the literature pUl CGI [13] are 
special cases of the above parametrization and satisfy these seven constraints. However, they 
are more constrained than the general parametrization obtained here. A special property of 
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and 



the TM mixing matrix of Eq. (11), also shared by its earlier parametrizations [10 
is that its first and third columns add up to zero. 

The first four of the seven unitarity constraints can be seen as normalization conditions and 
the last three as orthogonalization conditions. The three orthogonalization relations can be 
solved to get three real variables, say, -ReL^, ImU-zz and IvnXJ%\ which can be substituted 
back into the four normalization relations. This eliminates these three parameters. Two 
classes of solutions are obtained on solving the four normalisation relations in the remaining 
five parameters: one is CP conserving while the other is CP violating. 



5 The CP conserving solution 

The first solution has three independent parameters which can be chosen as ReUu, ReUis 
and ImUis. The imaginary part of Uu can be calculated from the equation 



|£/n| = V 2 /3-|^i 3 | 2 . (19) 

The elements U21 and C/23 are given by 

and 

The resulting mixing matrix has no Dirac type CP violation since first two elements in 
either first or third column have identical phases which results into a vanishing J. There 
are no constraints on the phases of Uu and U\$ and they are not necessary to describe the 
phenomenon of neutrino oscillations. However, the CP-odd rephasing invariants I\ and I2 
do depend upon these phases and are given by 

h = --ReUnImUu (22) 

and 

I 2 = Re{Ul x )Im{Ulz) - Im{Ul x )Re{Ulz). (23) 

This solution is called CP conserving as there is no Dirac type CP violation here. Tri- 
maximality is a mass-independent texture and it puts no constraints on neutrino masses 
and Majorana phases (Sj. Hence, a special case with vanishing Majorana phases will have 
real Uu and U^ parametrizable by a single parameter, say 9. The mixing matrix will be 
completely real in this special case and is given by 
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where c = cos 9, s = sin.9 and the parameter 9 can be related to the mixing angle #13 through 
the relation sin 6*13 = w| sin 9. This mixing matrix coincides with the mixing matrix obtained 
in a lepton mass model proposed by Grimus, Lavoura and Singraber [2] for the first sign 
choice. 



6 The CP violating solution 

The second solution has four independent parameters which can be chosen as ReUu, ReUi 3 , 
ReU2\ and ImUi 3 . The imaginary parts of Uu and C/ 2 i are calculated from the equations 



|C/ii| = V 2 /3-|^i 3 | 2 (25) 

and 

TmJIm = — 
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ImU 21 = -- {ImUn + D) (26) 



where the real parameter D is given by 

D 2 = 3\U 13 \ 2 - [Re(U n + 2U 21 )f (27) 

which can either be positive or negative. The element U 23 is given by 
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The Jarlskog invariant J is given by 



J = - (ReUnlmUn + 2ReU 2 iImU 11 + DReU u ) ■ (29) 

6 

The solution is called CP violating solution since J is non-zero here. The CP conserving 
solution can be recovered from this solution in the limiting case where J = 0. J vanishes for 
the two values of U23 given by Eq. (21) which are also the two limiting values within which 
U23 varies in the CP violating solution [Eq. (28)]. Thus, the CP conserving solution for U23 
defines the limiting region in which U23 varies in the CP violating solution. 

A notion of maximal CP violation naturally emerges here. It can be seen from the Eqs. (28) 
and (29) that J is maximal when 9 23 is maximal and the maximal value is given by 

</ 2 = ^ii| 2 |£A3| 2 . (30) 

It should be noted that the maximal CP violation does not mean that the CP violating phase 
5 will also be maximum when J and #23 are maximal. The maximal value of J obtained 
here coincides with the intrinsic CP violation in the mass model by He and Zee [12]. This 
is because the mixing matrix in the model is a special case of the most general unitary TM 
mixing matrix considered here. 



7 Numerical analysis and results 

The relationship between 6*12 and 6*13 given by Eq. (19) or (25) has been depicted in Fig. 1 and 
is identical for both the solutions. The relationship is well known and has been reproduced 
here to show the values of #12 and #13 which will enter into the numerical calculations through 
Uu and C/13. The value of #i 2 is constrained automatically to be well below its 3a range 
when #13 is varied upto its 3<r upper bound pp. 

Both the CP conserving and CP violating solutions have identical predictions for the ee 
element of the neutrino mass matrix 

M ee = miC/i 2 i + ^+m 3 [/i 2 3 (31) 

whose magnitude, also called effective Majorana mass, will be observable in neutrinoless 
double beta day experiments. Varying the lowest neutrino mass over a reasonable range 
(1CT 3 eV-1 eV), calculating other mass eigenvalues from the squared-mass differences Amf 2 
and Amf 2 (Table 1) sampled using normal distribution and using the values of #12 and #13 
sampled in Fig. 1, the predictions for \M ee \ have been shown in Fig. 2 for normal and 
inverted orderings of neutrino masses. This graph matches with usual plots of \M ee \ as the 
function of lowest neutrino mass and has been reproduced here merely to check correctness 
of the numerical analysis. 

The variation of If/23 1 with | C/13 ] for the CP conserving solution can be calculated directly 
from Eq. (21) as 

|£/ 2 3| = ^(V3|£/ii|±|£/i3|) (32) 

and has been plotted in Fig. 3. It is noted that It/23 1 depends only on \U\s\ for the CP 
conserving solution. However, U23 depends on ReUn, ReUiz, ImUis and ReUix in the CP 
violating solution. Varying these free parameters randomly in their allowed ranges, the 
variation of [ C/23 1 with |£/ 13 | for the CP violating solution can be calculated from Eq. (28) 
and has been shown in Fig. 4. One can easily see that the CP conserving solution depicted 
in Fig. 3 sets the upper and lower limits within which the CP violating solution shown in 
Fig. 4 can vary. The same fact is true for the variation between #23 and #13 depicted in Fig. 
5 for CP conserving solution and in Fig. 6 for CP violating solution. The two branches of 
CP conserving solution for #23 can be viewed as the bounds within which #23 varies in the 
CP violating solution. It can be easily seen that #23 can be maximal even for non-zero #13 
in trimaximal mixing. The allowed span of #23 around the maximal value increases with 
#13 . When #13 reaches its 3<r upper bound, #23 also spans its 3<r allowed range p. The span 
within which #23 is allowed to vary increases almost linearly with #13 showing a non-linear 
behaviour only at very large values of #13 which has not been depicted in Figures 5 and 
6. The deviation from maximal 2-3 mixing is zero when CP violation is maximal and it 
increases as CP violation decreases and, finally, there is no CP violation at the point where 
the deviation from maximality takes its maximal value allowed in TM mixing. 

For the CP violating solution, J, I\ and I<i are calculated from Eqs. (22), (23) and (29) 
by varying the four free parameters ReUn, -Ret/13, ImU-ys and ReU^i in their full ranges 
allowed by experimental as well as the unitarity constraints. The relationship between J 

8 



with If/23 1 for I C/13 1 = 0, 0.1 and 0.2 has been depicted in Fig. 7. The central point is the 
TBM solution for which both J and t/13 are zero and 1 1/23 1 — -75- However, there are values 
of It/23 1 other than the TBM value for which J is zero [Eq. (32)] and there are values of J 
other than zero for which It/23 1 takes its TBM value [Eq. (30)]. The relationship between 
Z/13, t/23 an d J shown in Fig. (7) is unique to trimaximal mixing not found in other mixing 
schemes deviating from TBM mixing j9] (and refernces therein) and subject to verification 
in future neutrino oscillation experiments. 

The relationship between 5 and #23 for #13 = 4, 8 and 12 degrees has been depicted in Fig. 
8. It is apparent from this figure that 5 does not attain its maximum value at #23 = |. So, 
the point of maximal CP violation is not the point at which S is maximum. The maximum 
of 5 occurs when #23 is slightly below the maximal value. Furthermore, the span in which 
#23 is allowed to vary scales up with 6*13 (see Fig. 6 as well) but the allowed range of 5 is 
almost unaffected by variation in 6*13 showing a marginal increase of about 10 degrees with 
increase in #13. The relationship between CP violation and deviation from maximality can 
be visualized in the most apparent way in Fig. 9 which plots S23 as a function of J for 
#13 = 4, 8 and 12 degrees. For a given #13, the allowed values of J and s^ form an ellipse 
around the TBM point (J = 0, s\% = ~). The semi-major axis of the ellipse is given by the 
maximal value of J [Eq. (30)] and the semi-minor axis of the ellipse is defined by the CP 
conserving values of S23 [Eq. (32)]. 

For completeness, the variation of J and I2 with t/13 has been depicted in Fig. 10. There 
are no constraints on I\. Since the phases of U\\ and t/13 are n °t constrained in both the 
solutions for TM mixing, there are no constraints on Majorana-type CP violating phases a 
and (3. 



8 Conclusion 

In conclusion, TM matrix has been defined as a unitary matrix in which the three elements 
of the middle column are equal and their phases can be rephased away by phase redefinition 
of lepton fields. This definition of strong trimaximality is not rephasing invariant and differs 
from the weak trimaximality which demands only the equality of magnitude of the elements 
of the middle column without putting any constraints on the phases and, consequently, is less 
restrictive. The unitarity constraints on the TM mixing matrix have been solved and two 
solutions have been obtained. The first solution is CP conserving with vanishing Jarlskog 
rephasing invariant J and has three free parameters viz. ReUn, -Ret/13 an d irn.t/13 which 
completely determine the mixing matrix. The deviation of #23 from maximality scales up 
with #13 and is within experimental constraints. The second solution, called CP violating 
solution, has -Ret/21 as one more free parameter and reduces to the CP conserving solution 
in the limiting case of vanishing J. The CP conserving solution defines the upper and lower 
bounds within which #23 varies with respect to #13 in the CP violating solution. So, J 
vanishes whenever # i3 is zero or #23 equals its upper or lower bound for the given value of 
#13. The Jarlskog invariant J takes its maximal value — ^?= 1 1/11 1 1 t/i3| when #23 is maximal. So, 
a notion of maximal CP violation naturally arises here: CP violation is maximal for maximal 



2-3 mixing. However, the Dirac-type CP violating phase S does not takes its maximum value 
at the point of the maximal CP violation. The maximal value of 5 occurs when #23 is slightly 
less than its maximal value. The allowed range within which S can vary shows a marginal 
variation of about 10 degrees with $13. For a given value of #13, s| 3 and J form an ellipse whose 
center is defined by the TBM mixing and whose axes are governed by the maximal values 
of CP violation and deviation from maximal 2-3 mixing. Trimaximality puts no constraints 
on the Majorana-type CP violating phases a and /3 since it is a mass-independent texture. 
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Figure 1: Relationship between 6*12 and 6*13. 
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Figure 2: The effective Majorana mass \M ee \ as the function of lightest neutrino mass mo 
for normal hierarchy (lower branch) and inverted hierarchy (upper branch). 
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Figure 5: Relationship of #23 with # 13 for CP conserving solution. 
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Figure 6: Relationship of #23 with #13 for CP violating solution. 



13 




-0.10 -0.05 0.00 0.05 0.10 
J 

Figure 7: Relationship between J and | C/23 1 f° r \U\z\ = 0,0.1 and 0.2. 




Figure 8: Relationship between S and #23 for #13 = 4, 8 and 12 degrees. 
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Figure 9: The variation of J with sl 3 for 6*13 = 4, 8 and 12 degrees. 
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Figure 10: The variation of J and I2 with | C/13 1 . 
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